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Abstract 

We construct a static axisymmetric wormhole from the gravitational field of two 
Schwarzschild particles which are kept in equilibrium by strings (ropes) extending 
to infinity. The wormhole is obtained by matching two three-dimensional timelike 
surfaces surrounding each of the particles and thus spacetime becomes non-simply 
connected. Although the matching will not be exact in general it is possible to make 
the error arbitrarily small by assuming that the distance between the particles is 
much larger than the radius of the wormhole mouths. Whenever the masses of the 
two wormhole mouths are different, causality violating effects will occur. 
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1 Introduction 



The field equations of Einsteins theory, being local, do not fix the global structure of space- 
time. Kurt Godel [||] noted first in 1949 that general relativity admits topologically non- 
trivial solutions exhibiting causal paradoxes, including closed timelike curves. Whatever 
one's attitude to these paradoxes - whether one views them simply as embarrassments, as 
a spur to re-examining the foundation of physics, as a challenge to philosophers of physics, 
or as a clue to new physical principles which might outlaw or tame the pathologies - the 
issues they raise are non-ignorable and have recently received a good deal of attention 
[2-5]. 

In 1988, M. Morris and K. Thorne conceived the idea that an advanced civilization 
might be able to construct a traversible wormhole which connects two distant regions of 
space. To prevent the throat from closing, exotic material (with negative energy-density) 
needs to be packed into the hole. Such a wormhole can be employed as a time-machine 
by setting the two mouths into high speed relative motion. V.Frolov and I.Novikov Q 
pointed out in 1990 that even a static wormhole functions as a time-machine if the two 
mouths are at different gravitational potentials, e.g. if one of the mouths is held near a 
neutron star. These are probably the simplest models for illustrating and testing ideas 
about causality-violating spacetimes. 

The two mouths appear as two separate masses in the asymptotically flat exterior 
space. Spacetime can therefore not be spherically symmetric. Moreover, the mouths 
will naturally gravitate towards each other, thus violating the assumed staticity. The 
nonstaticity can be made arbitarely small for sufficient large initial separation of the 
mouths. Alternatively, it can be eliminated altogether by anchoring the mouths, e.g. by 
strings held at infinity. It is the aim of this paper to develop an explicit, strictly static, 
wormhole model of this type. 

In Section II we start from the so-called Bach-Weyl solution which describes two static 
particles held in equilibrium by a strut. The strut may be replaced by strings that are 
attached to each of the particles and run off to infinity. In Sec. Ill we construct a worm- 
hole from the Bach-Weyl solution by cutting out the interior of two timelike tubes each 
surrounding one of the particles. Matching these boundaries gives rise to a layer of (non- 
standard) matter. The form of the boundaries is described by two " deformation- functions" 
which are required to make the induced metric continuous. We show that solutions for 
the deformation-functions exist (Appendix A) and give their explicit form for the case 
that the two wormhole mouths are far apart. However, this yields only two out of the 
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three components of the induced metric continuous. To smooth out the remaining metric 
component we introduce additional matter. The result is that this does not significally 
contribute to the stress-energy tensor as long as the mouths are sufficiently separated. In 
Sec. IV we calculate the stress-energy tensor for the surface layer. Finally, we point out 
that such wormholes may be used as a time-machine whenever the mass parameters of the 
two wormhole mouths are different. 

2 The Bach-Weyl-Solution 

In this section we introduce the so-called Bach-Weyl solution which describes two static 
Schwarzschild particles held in equilibrium by a strut between them. For this we follow the 
work of W. Israel and K.A. Khan Q who calculated the gravitational field of arbitrarily 
many collinear particles. Any static axially symmetric spacetime can be written in a 
cylindrical coordinate system /(T, r,z,ip) in the form 

depending on the two functions z^(r, z) and A(r, z). The vacuum Einstein equations require 
X{r, z) to be an axisymmetric Newtonian potential function of (r,z,(p), treated as cylin- 
drical coordinates on a fictitious flat background. It is known that if the Schwarzschild 
line element of a particle with mass m is transformed to the coordinate system / , A is 
formally the Newtonian potential of a uniform rod with length 2m. Therefore, we may 
use the following picture: Consider two non-overlapping rods placed along the z-axis, their 
centers located at z = Oj (i=l,2). We assume that the i^^ rod has mass 6, and define the 
quantities pf{r,z) which measure the distance from its ends (see Fig.l): 

{Pt? = r' + {ztf, {Pi? = r' + {zrf 
zf = z - {ai + bi), zl = z - {ai - hi) (2) 

Following Q we introduce the notation 

E{i+,j+) = ptpj + Z+Z+ + 
E{z-J+) = E{j+,r) = pTp+ + zrz+ + r2 (3) 

with an analogous expression defining E(i~ As a solution of the Einstein equations 
for our two body problem we choose A to be the Newtonian potential of these rods: 

A (r, z) = Ai (r, z) + A2 (r, z) , A, = i In t ^- ~ (4) 

2 PT + Pi +^bi 
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The corresponding z^(r, z) is given by the sum 

''^'^ A''' E{i+,j+)E{i-,j-) 



^ ^ 1 E(i~ j~^)E(i~^ j~) 

i=i j=i 



With this choice the constant of integration has been adjusted to make v vanish at spatial 
infinity. Thus the spacetime is asymptotically flat and satisfies the vacuum equations 
everywhere except on the segment of the z-axis between the the two rods. 

VQ = v[r = {),z) = —r\, 2; G ai + 6i,a2 - 62 (6) 

This phenomenon was interpreted by Bach and Weyl as a strut which holds the two bodies 
apart. In one of the authors has considered how a string energy density and tension 
can be defined for such line singularities. For our case the non-vanishing components of 
the energy-momentum tensor L^^ of the line source are given by: 

L[] = L2= l(i_e-o) >o if 61,62 >0 (7) 

The strut therefore has a negative energy density and a pressure numerically equal to it. 
Another possibility for keeping the two particles apart which avoids a strut with negative 
energy density is offered by a different choice for the constant of intergration in equation 
(^). If e.g. we choose v = —1^0, our solution would be regular on the segment of the z-axis 
between the two rods. On the other hand, there are now two singularities extending from 
the two rods to spatial infinity. Staticity requires that these "strings" have tension and 
positive energy density. 



3 Matching of Wormhole Mouths 

We cut out the interior of the surfaces Si and ^2 surrounding each of the two particles - 
the wormhole mouths - and match the surfaces to get a non-simply connected spacetime. 
Let us introduce two spherical polar coordinate systems / and //, one centered at z = oi 
the other at z = 02 (see Fig.l): 

I{t,ri,^,ip) : t = e^^T 

r = \J ri (ri — 26i ) sin {) 

z = ai + (ri — 61) cos 'd (8) 

r = \Jr2{r2 - 262) sin d 
z = a2 + {r2 - 62) cos (vr - i9) (9) 
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We do not distinguish time and angle coordinates of the different coordinate patches 
because we want to identify points on the wormhole mouths with equal values of (t, ??, ip). 
The constants Aj determine the transformation of the time coordinate of the metric (|l]) to 
the charts I and II. We will see that they cannot be equal in general which is the crucial 
fact leading to the occurence of closed timelike curves (CTCs). Notice that the polar axis 
of the spherical coordiate system near mouth 5*1 points to the +z - direction but the polar 
axis of the spherical coordinate system near mouth 5*2 points to the —z - direction. With 
this choice we identify the inner poles A and B of the wormhole mouths (see Fig.l). 
In these coordinates the Bach-Weyl-solution becomes 

_e2(A.-A.)(i_^)^,2^ .^^.^1^2 ; (10) 

where ct, is defined as the perturbation of the function i/(r, z) induced by the other particle, 
i.e. 

y = V22 + 2!/l2 = 1^11 + CTi = Z^22 + ^2 (11) 

Now let the parametric equations of the surfaces Si and 5*2 in our charts I and // be 
x-* = x-'(i^") with three-dimensional intrinsic coordinates = {t,'d,ip) - notice that i,j = 
1,2 but a, 6 = 0,2, 3: 

Si: ri= Re''-''^ 

S2: r2 = i?e^W (12) 

The deformation functions €{'d) and 6{'d) just depend on the angle and will be determined 
in the following. First we calculate the induced metric on Si and S2- The three holonomic 
basis vectors 6(^)1 5^ = -^\si tangent to Si have components ej^^^l^. = and their scalar 
products define the metric induced on Sii 



9ab\s, = e(a)|5, • e(fe)|5, = 9al3efa)e(^b)\Si (13) 



which for Si reads: 



To get the induced metric on 5*2 one has to interchange the indexes 1 and 2 and to replace 
5 by e. Next, we match these surfaces by identifying points with equal values of {t,i},ip), 
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i.e. Si = S2 = S. To obtain a continuous four-metric the induced metric has to be the 
same on Si and S'2. We may achieve continuity on S for the metric coefficients 500 and g^^ 
by specifying the functions 6{'d) and e(i?) and the relation between the asymptotic time 
scales, 

In general, continuity for 522 can not be achieved exactly. However, we will see that 
the error in the 522 component can be made arbitrarily small by enlarging the distance 
D = a2 — CLi between the wormhole mouths. 

Throughout this paper we assume that the mass parameters bi are positive. Hence, to 
make sure that the surfaces Si are timelike we impose the condition 

26i,2 <R^D (16) 

In Appendix A we show that under this condition also the radii Re^ and Re'' are greater 
than the Schwarzschild-radii 2bi and 262 respectively. 

To calculate the deformation functions 6 and e we equate the metric-coefficients g^slsi 2 
and Qoolsi 2 which yields a system of two equations: 

-Xiie,^) + e{^) = -X2{6,^)+6{^) (17) 



(,_|L).e^(,_||) as) 

There are different possibilties for fixing the size of the wormhole mouths, e.g. by choosing 
5 at an arbitrary point on the surface Si or by choosing the constant C. Notice that if the 
masses hi and 62 are equal, exact matching of the two wormhole mouths is possible. In this 
case it follows from equations ( [T^ and ( p!8| ) that 5(-d) = e{^) and C = 1. Such a spacetime 
will be symmetric with respect to the plane z = ^(02 — ai) and will not contain a Cauchy 
horizon. The identification of the mouths can be chosen in such a way that either CTCs 
exist throughout the whole spacetime - an eternal time machine - or that there do not 
exist CTCs at all. Under the condition (|l6| ) we can show that if bi 7^ 62 then C 7^ 1 in 
order that a solution to the system of equations ( p!7| ) and ( p^ ) exists. Hence, although our 
constructed spacetime is static, there does not exist a global timelike Killing vector field. 
In the terminology of V.Frolov and I.Novikov |p the gravitational field is nonpotential 
and according to their general proof CTCs have to occur. V.Frolov and I.Novikov studied 
explicitly a spherically symmetric wormhole model with two asymptotic regions. In order 
to be able to compare our results with theirs we choose C to have the value 



C 



(19) 
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With this choice of the constant C and under the assumptions (16) we are able to show that 
a solution 6 and e to the system (|l7|l8) exists for every ??. To calculate the deformations 
we eliminate (A2 — Ai) from eqn. (Il8|) by using eqn. (O) and solve for 6: 



(20) 



R \ \RJ 



I +C2|e 



2e 



262 

-r' 



Putting this result back into eqn. (|T^ leads to an implicit equation for £(■(?) which cannot 
be solved exactly: 

Fui^, e') := e' - e^i(^)-^2(5)g5(e) ^ q (21) 



However, if he wormhole mouths are infinitly apart i.e. D to infinity eq.(21) reduces to 

Foo{'&,e') -e^^'^ =0 (22) 



By virtue of equations (19) and (|20| ) one checks that e = is a solution. This is to 
be expected since deviations from spherical symmetry should tend to zero as the distance 
increases. In Appendix A we show that for sufficiently large D a solution to eq.(]2l|) always 
exists. Moreover it is shown that the solution is well approximated by Newton's method 
starting with e = as zeroth approximation. To first order in R/D the solution is: 



R 
D 



(1 



2^w _ 
R R' 



(23) 



An analogous result holds for the function e'^^''^ with the mass-parameters bi and 62 inter- 
changed. We see that the wormhole mouths are slightly deformed spheres of radius R and 
therefore our spacetime may be directly compared to the wormhole spacetimes treated in 



Now we focus attention on the metric component 522- We have to check the behaviour 
of the quantities ^ and oi under our assumption (p!6[). Using the implicit function 
theorem and expanding in powers of ^ (see Appendix A) we find 



M'd'd^ Ids 



(24) 



The discontinuity of the induced metric component (722 (see eqn.(|l4D) depends on the 
angle 'd and is given by: 



where l\g22{'&) 



522 152 - 522151 

,2(<72W-aiW)| 



^2(aiW-A2W)^2g25^^^^(^) 



( d^\2 ( dS_\2 



(I - ^) 



+ !}-{ 



(1 



261 



1} 



(25) 
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Using the definition of the functions (Tj(r, z) given by (11), inserting the functions Vij from 



equation (|5|) and expanding this in powers of R/D yields the approximation: 

A,..W~ -''""-3'"""'" ».n^tf + o(g) (26) 

Thus, for given bi/R and 62 /-R the discontinuity can be made arbitrarily small by choosing 
the distance D to be suitably large without changing the constant C. Causality violating 
effects will be unaffected by the distance of the wormhole mouths. 

The price for constructing a transversable wormhole is of course that non-classical 
matter must be present. In our case there is an infinitely thin shell concentrated at the 
surface S. Because of the discontinuity of the metric on this surface the Einstein-tensor 
is not well defined there, and in principle we cannot apply the thin shell formalism to 
calculate the stress-energy tensor. By a small modification of the metric - which physically 
means that we introduce some additional matter on one side of the shell - we circumvent 
this problem without strongly infiuencing the results obtained by following the standard 
shell formalism. We add a perturbation term to the component g^^ of the Bach-Weyl 
metric ( |To| ) in chart / which should exactly cancel the discontinuity of the induced metric 
on the surface S: 

gMn,^) - 62(^^1 (^^•'')-^^(^i'''))rf{l + Ag22{^)fa{ri - Re'^''^)} (27) 

where we have introduced a profile function fa{x). This function should be at least twice 
differentiable and satisfy the properties fa{^) = 1 and fa{x) = for x > a. For example 
we may take it to be of the form: 



fa{x) = e("^-^^) for < a 

= for !x| > a (28) 

The stresses and the energy associated with the additional matter surrounding the shell can 
be analysed by calculating the Einstein-tensor of the new perturbed metric. As expected, 
it turns out that the energy- momentum tensor is of the same order as the perturbation 
function, i.e 0{R^ / D^), but also depends on the first and second derivatives of the profile 
function which are proportional to 1/a and l/a^ respectively. Therefore the matter dis- 
tribution has to be spread over a region thicker than R/D to keep its density low. Notice 
that in general this matter will also violate the energy conditions. 
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4 Stress - Energy of the Wormhole Shell 

The surface stress-energy tensor Sab of the layer is linked to the jump [Kab] '■= K^ah^ ~^ab 
of normal extrinsic curvature across S K^^^^ are the extrinsic curvatures corresponding 
to the different imbeddings of S in charts I and II, each defined as 

Kab-=-n-Ve^^)e(^a) (29) 

The vector n is the normal to the surface S pointing from chart / to // and is normalized 
to one, e(a) are the three holonomic basis vectors of the surface S defined above. For 
non-lightlike surfaces the following distributional equivalent of Einstein's field equations 
holds: 

- 87r(5,, - \gabS) = [K^b] = K^'^ - i^f? (30) 

Notice that the stress-energy-tensor does not depend on the sign of the normal vector n 
because if it would point to the opposite direction, the two terms on the right hand side 
of equation (30) would have to be interchanged. For the unperturbed metric (p^) the 



components of n with respect to the different charts are: 



— -t- I 



n 



As expected the calculation yields that the energy density of the wormhole mouth mea- 
sured by any observer will be negative, pressures will be positive: 



^2 = Si^^\ -r== + -^r== I +0(^) (33) 




If we apply the thin shell formalism to the perturbed spacetime (27) the result for the 
stress-energy tensor of the shell does not differ from the calculation using the original 
metric up to order 0(i?^/D^). This results from the fact that the perturbation in the 
metric and its first derivative with respect to are of this order. (For the profile function 
( p8| ) the first derivative with respect to the ri vanishes on S.) 
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5 Conclusion 



V. Prolov and I. Novikov Q considered the general situation of static wormhole space- 
times where the gravitational field is non potential. They proved that the time gap for 
clock synchronization in the external space with respect to synchronization through the 
wormhole handle grows with time. The wormhole becomes a time machine as soon as this 
time gap is larger than the time needed for light propagation between the two mouths in 
the external space. In our case the creation of a time machine can be seen explicitly by 
studying the coordinate transformations 9). Notice that these transformations are cho- 
sen in such a way that the resulting spacetime is symmetric with respect to the spacelike 
hypersurface T = t = 0. We see that any other slice of constant external time T which 
enters mouth Si will reemerge at mouth 5*2 at a different time T' = CT. Hence, there 
exists a time gap ATgap between the wormhole mouths which is changing linearly in T: 

ATgap = \T' -T\ = \{C -1)T\ (34) 

For the case that the wormhole mouths are kept in equilibrium by strings extending to 
infinity (independent of the choice of the original (10) or the perturbed metric (27)) the 
external time interval AT/vg needed by a null geodesic propagating from the inner pole 
A of wormhole mouth Si to the inner pole B of mouth 5*2 is given by: 



ATng = D- i?(e^(°) + e^(°)) + 62 + &i 

D-fie'^W +62 

Re^(o) _ 2hi ] 




+2h 

\L> -Oi- 132 J ~ \ 

, ^ -^z , n ^ D-Re<^^ + bi-b2 \ 

Because this is the shortest path between the two wormhole mouths, CTCs exist if the 
time gap ATgap is larger than this time intervall ATng- By symmetry with respect to 
the hypersurface T = t = the constructed spacetime consists of two regions containing 
CTCs separated by a region without CTCs. If for example 62 > 61 and therefore C > 1, 
a causality horizon will appear as soon as the null geodesic running along the z — axis 
from point yl to i? is closed. This happens at the value of exterior time Tch, where Tch 
is given by 

ATng = ATgap =: {C - 1)Tch. (36) 
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Appendix A 



We prove that for large D there exists a deformation function e(i?) which satisfies equation 

Fd{^, e') := e' - e^^(''^^-^^(^'^^e^ = (37) 

Notice that is a function of e*^ given by eqn. (pO|) and that Ai is to be taken on ^2 while 
A2 on Si: 



,Ai(e,i?)-A2(5,i?) 



46i 



1 + 



462 



\ V Pt + Pi +2h J \ p+ + - 262 



(38) 



pf\s^ = {Re'{Re' - 262) sm{'&f + {D - [Re' - 62) cos(i?) T hlf}^ 

p±|s^ = {Re^{Re^ - 2bi) sm{^f + {D - [Re^ - bi) cos(i?) ± 62)^}^ (39) 

We show that e^^''^ can be approximated by using Newton's method and additionally that 
its derivative ^e'^^'''* exists and can be calculated by implicit differentiation of ^^('i?, e*^). 

Since for hi = 62 the function Fo{'d,e'') vanishes identically, we assume - without loss 
of generality - that b2 > bi > 0. Under this condition ^^(■i?, e*^) satisfies the following 
properties which tell us that for every fixed there exists at least one zero on the interval 

) = ^-e^-^^f >0 (40) 

Ffl(t?,l) = l-e^i-^2 <o (41) 

These two inequalities can easily be verified for sufficient large D by expanding e^^^^^ in 
powers of ^, 

By estimating the remainder of the expansion we actually can show that properties (^) 
and (|1|) hold for at least D > AR. 
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For large D the function F£){'d,e'^) is strictly monotonically decreasing and concave 
with respect to e^. This can be seen by considering F£){'d,e^) and its first and second 
derivative with respect to e*^ as a sequence of functions labeled by D. For every fixed -i? 
this sequence converges uniformly on / to a limit function which can easily be estimated: 

lim -^Fn{^,e^) = 1-C^ ^""l ~P < -hz^ ^ on I (43) 

(f 

\im -—.FD{^,e') <Q on I (44) 

A general theorem about the convergence of Newton's method (see for example Heuser I, 
p. 408) says that if conditions (40,41,4^ and ^) hold Newton's sequence, defined by 

converges to a unique zero Q'va. I when started from Cq = 1. We perform one iteration and 
estimate the error by 

^ ^1 < L^M^l^ = Ci f-V where := min |-^Ffl(i?, e^)| (46) 

The existence of such a constant Ci for large D follows from the facts that /x is bounded 
away from zero (see eqn.(|4^)) for D tending to infinity and that the first two terms of the 
expansion of the function i<£)(T?, ef) in powers of ^ vanish. Without exactly determining 
Ci we may conclude that e\ (eqn. ([2^) ) gives the zero C, correctly to order ^. 

It is easy to see that the partial derivative £^Fo{^,e'^) exists and on the other hand 
we have shown that for sufficient large D the partial derivative ^i<£)('!9, e*^) exists and is 
nonzero (see eqn.(^)). Therefore, it follows from the implicit function theorem that e'^^'^^ 
is differentiable and ^e.''^^^ is given by 

Expanding this quantity in powers of ^ leads to the estimate given by equation ([2^). 
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Figure 




Figure 1: The relation between the cylindrical coordinate system 7(T, r, z, if) and the two spher- 
ical coordinate systems I{t,r\,'d,ip) and II{t,r2,'&,<fi) centered at the wormhole mouths Si and 
^2 respectively is shown. The time and angle coordinates (T, t, ip) are suppressed and we have not 
distinguished angle coordinates (t?, ^p) of the different patches 7, I and //. The figure also makes 
clear the geometrical meaning of the quantities pf defined in Section 2. 
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